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Abstract. Particles are accelerated in cosmic sites probably under conditions very different from those 
at terrestrial particle accelerator laboratories. Nevertheless, specific experiments which explore plasma 
conditions and stimulate particle acceleration carry significant potential to illuminate some aspects of 
the cosmic particle acceleration process. Here we summarize our current understanding of cosmic particle 
acceleration, as derived from observations of the properties of cosmic ray particles, and through astronom- 
ical signatures caused by these near their sources or throughout their journey in interstellar space. We 
discuss the candidate-source object variety, and what has been learned about their particle-acceleration 
characteristics. We conclude identifying open issues as they are discussed among astrophysicists. - The 
cosmic-ray differential intensity spectrum across energies from lO^^eV to lO^^eV reveals a rather smooth 
power-law spectrum. Two kinks occur at the knee (:^10"'^^eV) and at the ankle {c^3 lO^^eV). It is unclear 
if these kinks are related to boundaries between different dominating sources, or rather related to charac- 
teristics of cosmic-ray propagation. Currently we believe that Galactic sources dominate up to lO^^eV or 
even above, and the extragalactic origin of cosmic rays at highest energies merges rather smoothly with 
Galactic contributions throughout the 10^^-lO^^eV range. Pulsars and supernova remnants are among the 
prime candidates for Galactic cosmic-ray production, while nuclei of active galaxies are considered best 
candidates to produce ultrahigh-energy cosmic rays of extragalactic origin. The acceleration processes are 
probably related to shocks formed when matter is ejected into surrounding space from energetic sources 
such as supernova explosions or matter accreting onto black holes. Details of shock acceleration are com- 
plex, as relativistic particles modify the structure of the shock, and simple approximations or perturbation 
calculations are unsatisfactory. This is where laboratory plasma experiments are expected to contribute, 
to enlighten the non-linear processes which occur under such conditions. 

PACS. 96.40 cosmic rays - 95.30 astrophysical plasma 



1 Introduction 

Particle acceleration in laboratory environments mostly 
is based on generating an electric field to do the nec- 
essary acceleration work on charged particles. In cosmic 
sites, however, plasma currents would quickly short-circuit 
any eventually-present electrical fields. A magnetohydro- 
dynamic view of cosmic plasma implies that no local elec- 
trical fields exist in the co-moving plasma reference frame, 
so no acceleration is obtained. How then is the strikingly- 
efficient acceleration of cosmic ray particles actually pro- 
duced? Obviously, cosmic sites can provide environments 
for particle acceleration up to energies above lO^^eV [Ij 
(see FigjlJ. This is more than three orders of magnitude 
more than what we can achieve in terrestrial accelerator 
laboratories: The LHC at CERN will accelerate protons to 
energies up to 12 TeV in its final stage. A cosmic ray par- 
ticle with an energy around lO^^eV (100 EeV) had been 
detected as early as 1962 [2], more events were then seen 
above 10,000 TeV (=10 EeV) where fluxes are as low as 



one cosmic ray per square kilometer and year. As a field 
of astrophysics, research of relativist ic-particle accelera- 
tion in cosmic sites had been initiated by Victor Hess, Al- 
bert Gockel, and Werner Kolhorster[3 : With their brave 
balloon measurements of ionizing radiation and its depen- 
dence with altitude in 1910-1913 up to 9 km altitude they 
demonstrated the existence of ionizing radiation impact- 
ing the Earth from outside. Combining and comparing 
astronomical findings about cosmic rays with particle ac- 
celeration studies in terrestrial laboratories is an essential 
part of our path to an understanding of the sources of 
cosmic rays. The extremely-high energies accessible only 
through cosmic rays provide a unique extension of the 
energy range being studied, albeit incurring the need for 
astrophysical modeling of the 'experimental setup'. 

It is the purpose of this paper to present our cur- 
rent understanding of cosmic particle acceleration and the 
nature of its sources, and to summarize the open issues 
herein, in order to cross-fertilize discussions among physi- 
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Fig. 1. The intensity spectrum of cosmic rays extends over 
more than ten orders of magnitude in energy with a surpris- 
ingly smooth spectrum I{E) ^ ^-2-7... 3.0 (^^^japted from 4 

cists on how particles can obtain relativistic energies, from 
the viewpoint of observational/experimental astrophysics. 



2 Cosmic Rays 

2.1 Measuring Cosmic-Rays 

After the recognition of cosmic radiation around 1910- 
1913, many experiments were set up to measure more de- 
tails about cosmic rays in the upper atmosphere of the 
Earth and in near-Earth interplanetary space. From subse- 
quent work of Robert Millikan, Arthur H. Compton, Pierre 
Auger, and others it became clear the ionizing radiation 
of extraterrestrial origins identified by the balloon experi- 
ments of Hess and others [3 [5 [6 [7 were caused by charged 
particles of extremely high energy and of possibly extra- 
galactic origin, which created showers of secondary parti- 
cles in the atmosphere. The discoveries of new elementary 
particles in cosmic ray showers, such as the muon, but also 
nucleonic excitation states, marks the decades of fruit- 
ful stimulation of nuclear physics and high-energy exper- 
iments by measurements on cosmic rays. Balloon experi- 
ments pioneered various techniques to record cosmic rays 
more directly than through their secondary showers, using 
emulsions and direct spectrometric analyzers with track- 
ing chambers and magnets. Nowadays space experiments 
such as STARDUST [8^ and AMS [9 implement the tech- 
nological forefront of those same technologies for captur- 
ing cosmic ray particles directly in interplanetary space. 
The space missions of the Advanced Composition Explorer 



ACE (launched 1997) [10^ and PAMELA [IT] (launched 
2006) provide detailed cosmic-ray intensity spectra dis- 
criminating different nuclear and lepton species (e.g. [12 ), 
and balloon experiments such as HEAT (1994) and ATIC 
(2000, 2003) have a role specifically in searches for anti- 
matter in cosmic rays |T3] [M] [15] . 

Above energies of lO^^eV, particle fluxes become too 
low for such direct observations. Demands on detector vol- 
ume and analyzing-magnet strength become excessive at 
energies exceeding thousands of GeV per nucleon. The 
technique of cosmic-ray studies then relies on using the 
Earth atmosphere directly as a detection volume, infer- 
ring the cosmic-ray particle characteristics from detailed 
measurements on the secondaries. Electromagnetic show- 
ers are initiated from the first interaction of an incident 
cosmic-ray particle at an altitude of 15 km and above. 
A variety of signatures can be recorded to re-constitute 
the nature (and direction) of the primary particle. Some 
of these techniques are only indirectly related to the pri- 
mary particle, such as the widths or leptonic contents of 
showers, some are more direct such as Cerenkov radia- 
tion caused by particle motion at a velocity in excess of 
the local- medium's velocity of light. Nitrogen fluorescence 
light is emitted by air along trajectories of shower secon- 
daries, and can be photographed from different viewing 
angles to construct a 3D shower image. The energy and 
flux of cosmic rays have been measured quite successfully 
with such techniques up to the extremes. Exciting results 
have been reported e.g. from the KASKADE ground level 
shower detector array of compositional changes around 
the knee of the cosmic ray spectrum [16 . But the nature 
of cosmic-ray particles is difficult to disentangle from the 
characteristics of the secondaries, as the first interaction 
with nuclei of the upper Earth atmosphere (or even inter- 
planetary gas) breaks up nuclei into the individual nucle- 
ons and erases the primary composition. Extensive Monte- 
Carlo calculations of shower properties are being used (e.g. 
[17]), but we must bear in mind that those calculations 
extrapolate high-energy physics interaction processes and 
cross sections significantly beyond where calibrations are 
feasible, and may not be realistic. 

A combination of techniques has proven most success- 
ful for the current generation of cosmic-ray telescopes at 
highest energies: Secondary shower particles are recorded 
in surface detectors, and several telescopes capture an im- 
age of fluorescence light emitted along the trajectories of 
particles through the atmosphere. The AUGER experi- 
ment implements this technique over a large surface area 
of 3000 km^ near Mendoza in Argentina [18], and a cor- 
responding site AUGER-North is under development in 
Colorado, USA. 



2.2 Indirect Observations about Cosmic Rays 

Cosmic rays interact with matter and radiation along their 
trajectories through interstellar space. This produces char- 
acteristic astronomical signatures, which add unique in- 
sights into the nature of cosmic ray sources. 
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Secondary gamma-rays are readily observed at GeV 
energies, now also up to the TeV range, from interstellar 
gas, which generates a source of photons byh such cosmic- 
ray interactions through the processes of Bremsstrahlung 
and inverse- Compt on emission as well as hadronic interac- 
tions [19][20]. Early measurements with SAS-II and COS- 
B established this aspect of Galactic cosmic-ray studies, 

which continued with the Compton Observatory and presentlf^econti knee) near lO^^'^eV and a c^3.3 beyond, before 



2.3 Characteristics of Cosmic-Rays 

The observed cosmic ray spectrum (Fig. [T]) extends over 
ten orders of magnitude of energy with a remarkably smooth 
spectrum in the form of a powerlaw I{E) = a • E~^. The 
powerlaw index a is c^2.7 at energies up to lO^^eV (the 
knee), then steepens towards a c:^3.0 with a small kink 



flourish from data captured with the instruments of the 
AGILE and FERMI missions. We distinguish several spec- 
tral regimes of the connection between cosmic rays and 
gamma-rays, depending on the characteristics of the pro- 
ductions of secondaries [21 . For ultrahigh primary cos- 
mic rays in low-density surrounding gas with weak mag- 
netic fields, secondaries of the electromagnetic cascade ini- 
tiated by collision of a cosmic-ray particle with ambient 
plasma constituents are mainly produced in forward di- 
rections, and carry the same order of magnitude in energy 
as the primary cosmic rays. If the magnetic field is strong 
enough to significantly affect the motion of secondary elec- 
trons at the low-energy tail of the electromagnetic cas- 
cade, isotropic lepton clouds and thus a pair-production 
halo around the source would accompany the high-energy 
secondaries - this halo, however, has not been measured 
yet. For even stronger magnetic fields, synchrotron energy 
losses result in a complete loss of the primary cosmic- 
ray's direction and energy, and lead to characteristic syn- 
chrotron emission caused by the high-energy cosmic rays 
[22] . Observations of synchrotron emission extend from X- 
ray energies down into the radio regime, and are a valu- 
able diagnostic to establish the existence of cosmic rays 
in specific source regions, due to the excellent astronomi- 
cal precision of such observations. At energies below TeV, 
two major interactions dominate the secondary gamma- 
ray emission: The creation of pions (rest mass energy / 
threshold is 140 MeV in the collision rest frame) produces 
gamma-rays efficiently through the 27 decay of neutral 
pions TT^, and Bremsstrahlung is efficiently emitted from 
cosmic-ray collisions with ambient nuclei at energies below 
GeV. 

An additional process arises from the interaction of 
cosmic-ray electrons with ambient photons, up-scattering 
these in energy to contribute to the gamma-ray emissivity; 
this is the inverse Compton scattering, with a character- 
istic steep decline towards lower energies [23] [24] . Photon 
fields for such interactions are provided from the cosmic 
microwave background as well as from stars and interstel- 
lar dust, but often more importantly from the cosmic-ray 
sources themselves through accretion disks or bright cen- 
tral sources and jets. Although understanding the inverse- 
Compton emission implies knowledge or assumptions about 
these photon fields, the characteristic shape of the inverse- 
Compton emission spectrum provides an important dis- 
crimination of leptonic versus hadronic origins of the sec- 
ondary gamma-ray emission from cosmic rays. Inverse- 
Compton energy losses upon ambient photons are relevant 
for cosmic-ray electrons (and positrons), and limit the dis- 
tance which they can travel through interstellar space to 
~kpc, thus implying an origin within our Galaxy. 



its slope fiattens again near 3 lO^^eV to a c^3.0 at the an- 
kle. Beyond, data become sparse, but another steepening 
appears above :^ 5 • lO^^eV due to the Greisen-Zatsepin- 
Kuz^min cut-off {GZK) [25 [26 . This cut-off is expected 
as the cosmic-ray protons above this energy exceed the 
pion production threshold when interacting with cosmic- 
microwave photons. A tentative detection of the GZK cut- 
off had been reported from the HiRes experiment [27] . 
and the AUGER results at highest energies [28] are also 
consistent with this GZK cutoff, although the (uncertain) 
spatial distribution of cosmic ray sources and their inho- 
mogeneities within the nearest few hundred Mpc also play 
a role for the spectral shape at the high end. 

The composition of cosmic rays appears dominated by 
protons and atomic nuclei, electrons make up about 1-2% 
of the total cosmic-ray ffux at GeV energies, and even less 
above [29 . We must note that our knowledge about com- 
position is best only at energies up to GeV energies; only 
a few pivotal composition measurements exist at much 
higher energies (see above). Positrons are observed at the 
level of 10~^-10~^ of the total, and anti-protons are two 
orders of magnitude lower in abundances still, at 1 GeV, 
yet with a rising spectrum so that their abundance be- 
comes comparable to positrons at 10 GeV. Composition 
measurements at higher energies are a lively topic, re- 
ffected by the investments into balloon and space instru- 
mentation, most notably the ambitious AMS project on 
the International Space Station ISS planning to place a 
large super-conducting magnet spectrometer in space [30 . 
One of the main drivers here is the search for antimatter 
and for decay products of dark matter particles, both ex- 
pected to appear at energies of 100 GeV and beyond. Re- 
ports about evidences for excess in fluxes of positrons [SU 
[32] and electrons p. 5 at such high energies have re-kindled 
discussions about dark-matter and SUSY signatures in 
cosmic ray spectra. The composition of elements and iso- 
topes at sub- GeV energies are very detailed and precise, by 
contrast [12] [33] [34]. Surprisingly, elemental abundances 
in cosmic rays are quite close to solar abundances over- 
all. With closer inspection, remarkable systematic differ- 
ences appear, however: The light elements of charges 3-5 
(Li, Be, B) are much more abundant in cosmic rays than 
in the solar system, also intermediate-mass elements be- 
low the iron peak are somewhat enriched in cosmic rays. 
This is understood from spallation reactions of energetic 
cosmic ray nuclei upon their collisions with ambient in- 
terstellar Hydrogen along their trajectories in interstellar 
space. Therefore, the abundances of all secondary isotopes 
produced by spallation reactions are enriched in cosmic 
rays with respect to the standard abundance distribution 
[35] . The secondary-to-primary isotope abundances can be 
translated to an estimated average residence time of order 
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10'' y, using nuclear spallation cross sections, an average 
density for interstellar matter, and cosmic-ray velocities 
[36j . In particular several radioactive isotopes, which have 
been found in cosmic rays at low energies, most notably 
-'^^Be, ^^Al, and ^^Cl, determine such cosmic-ray propaga- 
tion ages to 15±1.6 My [37]. 

Those inferences, however, rely on a rather simple prop- 
agation model for cosmic-rays, the leaky-box model. In this 
model, cosmic rays are accelerated at their (unknown) 
sources, and diffuse through interstellar space as mediated 
by turbulent interstellar magnetic fields. At high energies, 
leakage of cosmic rays from the plane of the Galaxy into 
the halo is expected as magnetic rigidity increases. Al- 
though the leaky-box model is quite successful for illus- 
trating several cosmic-ray properties such as the produc- 
tion of secondaries, cosmic-ray propagation is a complex 
superposition of many processes and boundary conditions, 
addressed best in numerical models (e.g. GALPR0P,[19] ^ 
Current modeling accounts for 3-dimensional distributions 
of cosmic-ray sources, radiation fields, and interstellar gas. 
The interstellar medium of the Galaxy with its complex 
and rapidly-evolving morphology (see below) may plausi- 
bly lead to local and episodic deviations from such time- 
averaged modeling (e.g. [38 ). 

At higher energies, species abundances can only be de- 
termined more indirectly (see above). The KASKADE ex- 
periment \W showed that He nuclei are the most abun- 
dant species near lO^^eV, and that rigidity cutoffs from 
interstellar magnetic fields suggest a dominance of heav- 
ier nuclei in the regime above the knee, in the transition 
region from Galactic origin of cosmic rays towards an ex- 
tragalactic cosmic-ray component, which probably dom- 
inates above energies of lO^^eV. Results are somewhat 
confusing, as they suggest on one hand that the effec- 
tive rigidity cutoff of cosmic rays around the knee (which 
is at lower energies for lighter nuclei) leads to a heavier 
composition with increasing energy, on the other hand the 
reported Silicon abundances do not show any cutoff signa- 
ture, and iron is only seen at highest energies and not at 
all at the lower energies where protons and helium domi- 
nate. Clearly, an assessment of systematic uncertainties is 
needed before conclusions should be drawn. 

At energies below 10 GeV, the cosmic-ray spectrum is 
seen to bend down from the powerlaw shape describing the 
higher energies. This results from shielding by the mag- 
netic fields of the Earth and Sun, deflecting cosmic rays 
due to their charge [40 . Additionally, secondary particles 
created by cosmic-ray interactions within the Earth atmo- 
sphere partly escape from the atmosphere, and then may 
be trapped in the Earth's magnetosphere. In particular, 
secondary neutrons streaming from the top of the atmo- 
sphere will decay into protons and abundantly be captured 
in radiation belts [41 . Solar- wind particles will provide an 
additional admixture to the low-energy spectrum of cos- 
mic rays observed near Earth, so that inference of the ex- 
trasolar cosmic-ray contributions below 10 GeV/nucleon 
are uncertain. Low-energy cosmic rays are trapped by the 
Earth's magnetic fleld and move along its field lines. Thus 
charged-particles penetrate deep into the atmosphere to- 
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Fig. 2. Particle acceleration in an interstellar shock region: 
The feedback of the relativistic particle acceleration leads to 
substantial modifications to the structure of the shock in the 
the ory of diffusive shock acceleration (DSA) (figure courtesy 
H)son, 2007) 



wards the Earth's magnetic poles. Arctic light is the obser- 
vational evidence of the abundance of low-energy cosmic 
rays captured in those radiation belts^ with cutoff rigidity 
depending on geomagnetic latitude. 

This sensitivity to magnetic fields explains why an 
astronomy with cosmic rays is impossible: the curvature 
of the trajectory of a cosmic ray, as estimated from the 
Larmor radius (typical interstellar magnetic fields have 
strengths of a few /iGauss), even at a cosmic-ray energy 
of a PeV would be smaller than the distance to even the 
most nearby stars. So, except for the highest energies, di- 
rect astronomical traces back to their original sources are 
not possible from observations of cosmic rays. 

Generally, cosmic rays are found to arrive near Earth 
with an isotropic distribution, at energies well above the 
modulation by magnetic fields of Earth and Sun. There 
had been some reports about deviations from isotropy, i.e. 
clustering of directions for highest energy cosmic rays from 
the AG AS A experiment (see discussion by [l^). More con- 
vincingly, for events above 40,000 TeV the AUGER ex- 
periment recently reported even an astronomical result: 
Their trajectories point back to nearby active galaxies 
(distributed along the supergalactic plane which describes 
matter distribution of the local universe on the scale of 
~100 Mpc), a result based on 81 cosmic-ray events \2S\ . 



3 Astrophysical Environments for Particle 
Acceleration 

As we have seen above, energetic particles interact very 
efficiently with even low-density gas of the Earth's upper 
atmosphere. Hence, the energy for acceleration plausibly 
is provided either by macroscopic streaming or turbulence 
of interstellar gas, or within the energetic jet outffows initi- 
ated by violent accretion or explosion processes connected 
to compact objects and black holes. Interactions of ener- 
getic charged particles with such streaming gas or plasma 
produce what we observe as cosmic rays. 
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Astrophysical plasma conditions are characteristically 
different from the situation of an ideal gas (which often 
is in our minds when we attempt to understand cosmic 
gas). A major difference is the predominance of Coulomb 
interactions with their long range, as opposed to short- 
range forces which characterize collisional equilibration 
{thermalization) in an ideal gas. This implies that cosmic 
plasma is inherently viscous^ and the motions of electrons 
and ions are coupled over considerable ranges through 
electromagnetic field morphologies and dynamics; mag- 
netic fields play a dominant role in mediating kinematic 
differences of cosmic particle populations. When collisions 
between individual charged particles even become negligi- 
ble with respect to the electromagnetic field interactions, 
the plasma is called collisionless - the standard case for 
astrophysical plasma [24 [43 . 

This leads to considerable complexities in treatments 
on smaller scales where energy transfers are substantial, 
such as in interstellar shocks (see below). On larger scales, 
it is often a useful approximation to consider the mag- 
netic fields as frozen (or rather coupled) to the moving 
plasma. An assumed energy equilibration leads to the con- 
cept of a multi-phase interstellar medium in pressure bal- 
ance^ specifically cosmic-ray pressure balancing magnetic 
and interstellar gas dynamic pressures. Typical pressure 
values are ~10~^^dyn cm~^ with a Galactic-plane mag- 
netic field value near 2-5 fiG. 

The more modern view of the interstellar medium is 
that of a dynamical and transitory medium, where energy 
injections by winds and supernova explosions are never 
balanced on scales below :^ 100 pc (e.g. [44 ). Therefore, 
in astrophysical plasma equilibria often are not obtained, 
hence perturbation approximations are inappropriate for 
treatment of dynamical processes such as acceleration. 
For example, the interactions between charged particles 
of very different masses, the electrons and ions, exchange 
very little energy, compared to interactions between parti- 
cles of similar mass. Hence, thermalization occurs on dif- 
ferent time scales for electrons and ions, and yet differ- 
ent between electrons and ions of a plasma. (For exam- 
ple, in a supernova remnant, electrons are characteristi- 
cally thermalized within 100 y, while ion thermalization 
takes 4000 y, and electron-ion equilibrization needs 2 10^ y 
[24J ; supernova remnants are expected to efficiently accel- 
erate particles while they are young and magnetic fields 
in the shock front are high due to the steepness of the 
shock, however.) This limits application of magnetohy- 
drodynamic theory (MHD), which is based on the overall 
fluid properties of a plasma interacting with internal and 
external magnetic fields; mixtures of three different fluids 
with their respective properties add considerable complex- 
ity to such MHD descriptions. Note also that MHD theory 
is non-relativistic. Far from equilibria, plasma dynamics 
exhibits turbulence [24 . In cosmic-source environments 
plausibly connected to relativistic acceleration, injection 
of plasma turbulence is expected, as one deals with cos- 
mic explosions and formations of jets. The description of 
turbulent processes does not allow direct kinematic solu- 
tion of plasma flows, it is largely restricted to following 



the energy transport between different spatial scales of 
turbulence, an equilibrium turbulence [2T. 

These considerations illustrate why processes of rela- 
tivistic particle acceleration are commonly treated in ap- 
proximative models. Acceleration in shocks thus was de- 
scribed in the test particle approach^ where the cosmic- 
ray particle was injected into the externally-defined shock 
region with its density, magnetic-field, and velocity prop- 
erties derived from hydrodynamical analogies (e.g. [45 ). 
Calculation of the Fermi acceleration process then involves 
the determination of the mean free path of the cosmic-ray 
particle between scatterings on magnetic- field turbulence 
at the two different sides of the shock with their different 
bulk velocities {momentum convection process). Per scat- 
tering between the different sides of the shock, the energy 
gain is approximately dE/E ^ 2dv/v (i.e., Fermi- 1 pro- 
cess, to first order in velocity difference across the shock; 
stochastic acceleration from random-motion scatterings 
{momentum diffusion, Fermi-II-process) would show a quadratic 
relation). The energy flow from the shock region to the cos- 
mic rays is ignored in this approach by definition. In this 
test particle model, the interaction zones which lead to ac- 
celeration through scattering of the test particle are set up 
independent of this acceleration process. Hydrodynamics 
leads us to expect temperature and density jumps across 
a shock front, which can be calculated analytically (e.g. 
density contrast by factor 4, from the Rankine-Hugoniot 
relation) [45 . 

Refined theories of shock acceleration have been devel- 
oped. Currently, the theory of diffusive shock acceleration 
(DSA) provides a description for supernova shocks, which 
reflects rather closely the observations of shock region ra- 
diation properties [46] [47] [48] . In DSA, the feedback of 
relativistic-particle scattering on the overall plasma dy- 
namics is included, leading to revised shock structure in 
temperature and density (e.g. [49 ). For example, the large, 
step-like temperature jump across the shock is reduced in 
magnitude, and smeared over a larger region on both sides 
of the shock. Conversely, the density changes across the 
shock become considerably stronger up to a factor eight. 
Since the magnetic-field strength is related to density, this 
incurs amplifications of the magnetic field which were ab- 
sent in the test-particle treatment, but seem to occur in 
supernova remnants such as Cas A (see below). Note that 
this diffusive shock acceleration model also emphasizes a 
difference to typical laboratory situations, where magnetic 
diffusion may occur, while in astrophysical situations the 
magnetic field is frozen into the plasma and its high con- 
ductivity transfers diffusion processes to plasma currents 
only (cmp. [24^). 

After this general discussion, we now turn to specific 
cosmic environments where relativistic-particle accelera- 
tion is either plausible or has been seen already. 



From observations combined with astrophysical theory work, 
several candidate sites for production of relativistic charged 
particles have been identified. Although one or few of these 
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Fig. 3. Detailed models have been calculated for the emission 
of different types of electromagnetic radiation from a young 
supernova remnant accelerating particles. For the 350 year-old 
SNR Cas A at 3.4 kpc distance, most detailed observational 
constraints are available, for tesing those model predictions. 
This broad-band spectrum shows the different emission pro- 
cesses discussed in the text: (from low to high energies) Syn- 
chrotron emission is observable in a broad band from radio to 
X-ray energies; Bremsstrahlung occurs at hard-X and gamma- 
ray energies in the MeV to GeV band, pion-decay emission 
occurs in the GeV to TeV band as a broad flat-topped hump, 
similarly towards TeV energies the inverse- Compt on emission 
has its peak, (adapted from [46]; see also [501151] ) 

source types may dominate, this is still an open field of 
research. 



4.1 Supernova Remnants (SNR) 

The standard argument for supernova explosions as origin 
of the cosmic-rays arises from energy considerations: Prop- 
agation characteristics of cosmic rays imply that at high- 
est energies they cannot be contained within the Galaxy, 
but escape. Observing a powerlaw spectrum and an as- 
sumed steady state demands that sources must replen- 
ish the energy lost from the Galaxy in the form of ener- 
getic particles. With an estimated confinement time (see 
above) of 10^ years, and a local cosmic-ray energy den- 
sity of 1.8 eV cm~^, one estimates an energy loss of 
about 2 lO^^erg s~^. This energy loss is large by astro- 
physical standards, since most conventional cosmic ob- 
jects have luminosities far below this. Supernovae occur 
at a rate 1/(30-50 years), and eject a total kinetic energy 
on the order 10^^ erg, hence generate a total energy input 
of lO^^erg s~^. Although an efficiency of 10% for the ac- 
celeration of particles appears rather high, in general this 
energy balance argument is adopted to make SNR the best 
candidate source of cosmic rays. 

The acceleration may plausibly occur in the shock re- 
sulting from the explosion into surrounding gas. The dis- 
continuity between the exploding envelope traveling at 
several 1000 km s~^ and ambient interstellar gas at rest 
sets up a step change in density and temperature near this 
outward-propagating discontinuity. Charged particles may 
now cross the discontinuity, and may repeatedly scatter 
on either side to cross the shock again and again, thereby 
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Fig. 4. Convincing evidence of relativistic-particle acceleration 
within a young supernova remnant has been accumulated for 
the SNR J1713. 7-3946: X-ray emission attributed to electron 
synchrotron emission (small color image) appears spatially cor- 
related with gamma-ray emission at TeV energies (larger im- 
age), the latter expected to arise either from inverse- Compt on 
emission of the same electron population on ambient photons 
emitted by the SNR, or by pion decay from cosmic-ray protons 
accelerated together with the electrons. The broad-band spec- 
trum in the right part of the Figure illustrates the different 
components of expected radiation (adapted from [52]) 

gaining net energy according to the Fermi-I acceleration 
mechanism. Such acceleration is only limited by increasing 
Larmor radii as the particle gains energy. Yet, acceleration 
depends on the efficiency of reflecting scatterings across 
the shock, hence on shock density, and on the velocity 
gradient across the shock. Therefore, supernova remnant 
acceleration of relativistic particles is intrinsically time- 
variable, depends on the age of the supernova remnant 
and its surrounding density structure. Young supernova 
remnants at ages 100-10000 years are best candidates, in 
general. 

The smoking gun signifying that particle acceleration 
occurrs in SNR was believed to be the observation of syn- 
chrotron emission from relativistic particles in the mag- 
netic field of the shock region, best coinciding with other 
clearly non-thermal radiation from high-energy particles 
such as Bremsstrahlung, pion-decay gamma-ray emission, 
or inverse-Compton emission [53]. SNR 1006 was a well- 
known X-ray source. When TeV gamma-rays were reported 
from the Cangaroo experiment [54], it was believed that 
this indeed was the smoking gun. Later more sensitive 
measurement by HESS demonstrated that the Cangaroo 
result was in error, however, and TeV emission would be 
much lower or even absent. Instead, SNR RXJ 1713.7- 
3946 is now considered a well-established proof of locally- 
accelerated particles within this young supernova remnant 
(see Figjij [5l][55l [56] . Here the morphology of the X-ray 
emission seen with the Japanese ASCA satellite instru- 
ment follows closely the morphology of TeV emission 
as clearly mapped by the HESS telescopes. Unclear re- 
mains the origin of the TeV emission: If magnetic-field am- 
plification within the shock should be large, synchrotron 
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emission would be stronger and the corresponding electron 
population would be insufficient to explain the inverse- 
Compton emission. In that case, hadronic origin through 
pion decays from proton interactions in the SNR gas would 
be the explanation. But as the magnetic fields are uncer- 
tain, both explanations are possible. 



4.2 Stellar Wind Interactions 

Most massive stars are known to occur in binaries, i.e., 
massive stars interact with other stars formed within the 
same parental molecular cloud, and capture collisions re- 
sult in binary systems. Subsequent stellar evolution will 
lead to phases where stars have strong winds with mass 
losses on the order of 10~^ M© y~^. For a subset of the 
binary systems, 0-star winds may coincide with Wolf- 
Rayet wind phases of their companions. Depending on 
the orbital separations of such binaries, interacting stel- 
lar winds with typical wind velocities up to and exceeding 
1000 km s~^ set up conditions very similar to what has 
been described above for supernova remnants. Therefore, 
such sources have also been discussed as candidate accel- 
erators of cosmic-ray particles (e.g. [58 ). 



4.3 Pulsars 

Pulsars are rapidly rotating neutron stars, with rotation 
periods as short as milliseconds for objects with a mass 
comparable to the mass of the Sun yet much more com- 
pact with typical radii of :^10 km. In general, the magnetic 
axis of the neutron star will not be aligned with the axis of 
rotation. The magnetic field of the pulsar will rotate with 
the star, but as it extends further out, an extreme con- 
dition arises at a radius where field lines corotating with 
the star would reach the speed of light. As a consequence, 
field lines not closing between the neutron star's magnetic 
poles within this light- cylinder will remain open. Plasma 
moving along these field lines will escape from the mag- 
netosphere. Such loss of charges incurs plasma currents as 
well as regions devoid of charged particles further in along 
those open field lines. These devoid regions are called gaps, 
and are assumed to constitute the electrostatic accelera- 
tion regions for relativistic particles ejected by the pulsar. 
Radiation effects such as pair-production cascade emission 
and curvature radiation will produce X- and gamma-rays, 
relativistic aberration will add complexity and may be re- 
sponsible for caustics (e.g. [59 . 

Aspect geometries may differ between different pul- 
sars. Therefore, the association of peaks in the light curves 
over a pulsar rotation period cannot easily be assigned 
to either different poles of the rotating neutron star or 
different gaps sweeping over the observer's viewing direc- 
tion. It is unclear where exactly these acceleration gaps 
appear: close to the neutron star surface would be the po- 
lar gap, closer to the light cylinder the outer gap, and slot 
gap models describe intermediate cases, where the gaps 
would be very extended and located near the last closed 




p~Mm5EC p~i50mSEC p - 89 mSEC p-iosmSeC P~3emSec P~ 337 mSEC P- IflT mSEC 
TIME IN FRACTIONS OF A PULSE PERIOD 



Fig. 5. Pulsar geometries are complex, the different aspect 
angles make it difficult to assemble a coherent and generalized 
pulsar emission model. Due to the rotation of the neutron star 
with its magnetic axis inclined, pair plasma generated along 
the inner field lines will flow out at open field line regions near 
the light cylinder. Voltage gaps are expected to result further 
in within the magnetosphere, which accelerate particles to rel- 
ativistic energies up to possibly lO^^eV. It is unclear where 
exactly these acceleration gaps appear; the sketches show the 
gap viewing angles and light curves over the pulsar rotation 
phase for (from left to right) the polar gap, the outer gap\ and 
the intermediate slot gap models, respectively. The lower part 
of the Figure shows observed light curves for seven pulsars, 
where emission is seen up to gamma-ray energies, (adapted 
from 60 61 ) 

field line. Fig. [5] shows those geometries, and the associ- 
ated light curves expected under those viewing and gap- 
location cases. For comparison, the lower part of FigjS] 
shows observed light curves for seven pulsars, from which 
detection of high-energy gamma-rays provided clearest ev- 
idence of currently-ongoing relativistic particle accelera- 
tion. 



4.4 Microquasars and Active Galaxies, and 
Gamma-Ray Bursts 

Accretion onto compact stars and black holes is a well- 
known energy source and driver of high-energy emission 
of binaries at X- and gamma-ray energies. For a subclass 
of such binaries, the ejection of plasma jets with relativis- 
tic velocities was concluded from radio images taken at 
different epochs which showed apparent motion of radio 
emission knots at velocities in excess of the speed of light 
[62 . More gravitational energy is available when the com- 
pact star is not a neutron star with a surface terminating 
inflow at a radius of 10-15 km, but rather a black hole. The 
Schwarzschild radius of a stellar object with the mass of 
the Sun is km, correspondingly smaller for more mas- 
sive stellar remnants. Note that massive stars are created 
with a power-law distribution in masses extending out to 
about 80-100 Gravitational energy generated outside 
the Schwarzschild radius may escape the black hole and 
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contribute to energy transfers taking part in the genera- 
tion of these bi-polar plasma jets. About 30% of X-ray bi- 
naries with radio emission were found to exhibit jets, and 
about eight micro-quasars are established to exist within 
our Galaxy. Clearly, ejection of plasma jets also consti- 
tutes a smoking gun of a relativistic-particle accelerator 

m- 

Angular-momentum transfer in the flow of material 
from near such a compact object onto its surface results in 
formation of an accretion disk. Plasma flow within such an 
accretion disk may set up a dynamo, and thus be the driver 
of plasma jets perpendicular to the disk. Similar physical 
processes are expected, although at different scales, from 
accretion processes on supermassive black holes in the cen- 
ters of active galactic nuclei [64 . Here, black- hole masses 
are on the order of 10^ M©, and accretion is assumed to oc- 
cur from nearby stars, being disrupted near the black hole. 
With correspondingly-larger Schwarzschild radii, tempo- 
ral variations of the accretion luminosity will be slower and 
more difficult to trace observationally. Therefore, micro- 
quasars are the objects where astrophysicists strive to ex- 
plore details of the physical processes at the heart of these 
accreting black holes. 

Gamma-ray bursts of the long class (i.e., with burst du- 
rations in excess of 2 s) have been associated with collap- 
sars, the gravitational collapses of massive stars to black 
holes at the end of their evolution [65 . The late stages 
of such a collapse suggest a similar situation, whereby 
the compact core of the massive star is accreted by the 
centrally-forming black hole. Again, formation of an ac- 
cretion disk is expected, though with a lifetime on the or- 
der of milliseconds only, but the physical processes could 
be quite similar to what has been discussed in the context 
of jets formed in microquasars and active galaxies. 

In theories of gamma-ray bursts, the internal-shock 
model had been proposed to explain the generation of 
gamma-ray emission from within the plasma jet [66 . This 
model is quite successful in explaining the properties of the 
prompt gamma-ray emission and its variability. Plasma 
ejection in the jet with Lorentz factors of several hundred 
shift the 7 — 7 pair production threshold to sufficiently- 
high energies, so that the accelerated-particle radiation 
may escape, as observed. 

In principle, all these jet sources may accelerate par- 
ticles to energies beyond the knee, thus being candidate 
sources for significant contributions to the observed cosmic- 
ray flux. AGN and GRB jet sources are powered by the 
largest energy reservoirs, and hence may be responsible 
for the ultrahigh cosmic rays which reach us from extra- 
galactic sources. 

5 Open Issues 

Astrophysical first-order models are successful to explain 
much of the radiation observed from above candidate sources 
of cosmic rays. But some ad-hoc assumptions and first- 
order approximations are unsatisfactory, in particular the 
assumption of equilibria and the application of the per- 
turbation method. Therefore, the connections derived be- 



tween the source models the overall production of cosmic 
rays through such propagation models are highly uncer- 
tain. Derivations in particular of cosmic-ray fluxes and 
energy spectra, and of specific-source contributions, are 
vague. 

Among the issues being discussed are: 

— What are the maximum- achievable energies for differ- 
ent candidate cosmic-ray accelerators? Is there a tran- 
sition between different source types? 

— What are the acceleration efficiencies for different can- 
didate accelerators? Which magnetic-field amplifica- 
tions occur in real shocks? 

— How can the injection problem be solved? Does self- 
injection happen? (Are non-thermal tails sufficient?) 

— How critical is the geometry of shock regions? Shocks 
are presently assumed to be planar. How does one 
properly treat non-planar shocks and their losses? 

— Are observed GeV-TeV gamma-ray emissions leptonic 
or hadronic in origin? What are the constraints on 
magnetic fields within the shock region, on densities, 
and on the radiation fields available for inverse-Compton 
scattering? 

Obviously, the issues related to shock properties are all 
in some form accessible to laboratory experiments with 
plasma accelerated through laser pulses. Particle- in- cell 
(PIC) codes for numerical treatments of the complex plasma 
processes with particle acceleration are now becoming more 
realistic with increased resolution in 3D, making use of 
advanced computing architectures [67J. Still, how labora- 
tory conditions can be set up to perform a sufficiently- 
controlled physics experiment, this remains a major chal- 
lenge. Discussions such as experienced at this conference 
are steps in this direction. 
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